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Numerous experimental and theoretical studies have established that intramolecular vibrational 
energy redistribution (IVR) in isolated molecules has a hierarchical tier structure. The tier structure 
implies strong correlations between the energy level motions of a quantum system and its intensity- 
weighted spectrum. A measure, which explicitly accounts for this correlation, was first introduced 
by one of us as a sensitive probe of phase space localization. It correlates eigenlevel velocities with 
the overlap intensities between the eigenstates and some localized state of interest. A semiclassical 
theory for the correlation is developed for systems that are classically integrable and complements 
earlier work focusing exclusively on the chaotic case. Application to a model two dimensional 
effective spectroscopic Hamiltonian shows that the correlation measure can provide information 
about the terms in the molecular Hamiltonian which play an important role in an energy range 
of interest and the character of the dynamics. Moreover, the correlation function is capable of 
highlighting relevant phase space structures including the local resonance features associated with 
a specific bright state. In addition to being ideally suited for multidimensional systems with a 
large density of states, the measure can also be used to gain insights into phase space transport 
and localization. It is argued that the overlap intensity-level velocity correlation function provides 
a novel way of studying vibrational energy redistribution in isolated molecules. The correlation 
function is ideally suited to analyzing the parametric spectra of molecules in external fields. 



I. INTRODUCTION 

The nature of the dynamics of a molecule in a highly 
excited rovibrational state has been a subject of consid- 
erable study and debate for a number of yearsEJ. De- 
ciding the fate of a localized excitation in a molecule 
at high energies in terms of the time scales, pathways, 
and final destinaikaiJ|iaS||POsed a significant challenge for 
experimentalistsBMOT and theoristsBfeiy alike. In- 
tramolecular vibrational energy redistribution (IVR), as 
it has been called, is perhaps one of the most important 
phenomenon in chemical physics. It is now well estab- 
lished that IVR plays a crucial role in our understanding 
of reaction rate theories and the ability to effect mode 
specific chemistry. The ad-veni; of modern high resolution 
spectroscopic techniquestimD'El has led to a fairly detailed 
study of a number of molecules in highly excited energy 
regions. Such beautiful and precise studies have demon- 
strated the richness and complexity of IVR in various 
molecules and are an essential step towards realizing the 
goal of quantum control of chemical reactions. 

As a result of a number of experimental and theoreti- 
cal studies, a stani^djjasadiam for IVR called the tier 
model has evolvealBlja'B&QQla. In this model an opti- 
cally bright state is coupled to a bath of dark states in 
a hierarchical fashion. In a typical experimental situ- 
ation, the bright state corresponds to the eigenstate of 
the zeroth-order molecular Hamiltonian Hq which car- 
ries almost all of the oscillator strength. The form for 
this zeroth-order Hamiltonian involves the rigid rotor- 
harmonic oscillator part corrected perturbatively for non- 
resonant anharmonic interactions. The rest of the eigen- 
states of Hq having negligible oscillator strength are 



termed dark states. At low energies the excitation is 
localized since the bright state is not coupled or very 
weakly coupled to the dark states. However, at higher 
energies the bright state is coupled to the dark states 
by various perturbations Vj which include anharmonic 
resonances, coriolis and centrifugal couplings. These off- 
diagonal couplings are responsible for the flow of en- 
ergy from the initial bright state into the various dark 
states. The system is more appropriately described by 
the Hamiltonian Hq + Y]j Tj Vj in such energy regimes 
with the Tj denoting the strengths of the perturbations 
Vj. Thus, the experimental signature of IVR in the fre- 
quency domain is revealed as a fragmentation of a single 
rovibrational transition into a multiplet. Ample evidence 
exists in the literature indicating that the bright state is 
not coupled democratically to all of the possible dark 
states. Instead there is a small subset of the dark states, 
termed as belonging to the first tier, that strongly couple 
to the bright state. This first tier of states controls the 
initial stages of the energy flow and hence corresponds 
to the shortest timescale in the system dynamics. They 
are coupled to another set of dark states forming the 
second tier and giving rise to a second, longer timescale 
for energy flow. This pattern of the couplings leads to 
a hierarchical tier-model description of IVR in isolated 
molecules-. For instance, Callegari et. al. have recently 
observedE3 seven different timescales ranging from 100 fs 
to 2 ns for IVR from the first CH strech overtone of the 
benzene molecule. 

One of the key issues in IVR studies is the nature of 
the states that make up the various tiers. Identifying the 
nature of the dark states is an important first step to- 
wards the ultimate goal of determining the tier by tier 
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classification of the states, i.e. deciding the potential 
IVR pathways from the bright state of interest. Theo- 
retically, in order to know the destination of the initial 
excitation, it is necessary to identify the terms in the 
molecular Hamiltonian which are responsible for the cou- 
pling between the bright and the dark states. A molec- 
ular Hamiltonian is needed, but it is well known that 
determining the ab initio potential surface to sufficient 
accuracy is not routine except for small molecules. In 
order to make progress, a slightly different viewpoint has 
been adopted in the literature. Starting from the exper- 
imentally measured levels spanning the energy range of 
interest, an effective spectroscopic Hamiltonian is deter- 
mined which reproduces the levels to a certain accuracy. 
The anticipation is that such an effective Hamiltonian 
can be used to study the eigenstates and dynamics of 
the system in the given energy range. Effe^pB-spectr£fc. 
scopic Hamiltonians have been determinedt§Latpli30'ta 
for a number of molecules like phosphaathyneEj (HCB-L 
hypocholorous acidliJ (HOCl), acetyleneO, and waterEil 
to quote a few. The effective Hamiltonians typically in- 
volve anharmonic resonances which couple the zeroth- 
order modes of Hq. Fairly detailed quantum, classical, 
and semiclassical stud i es , u s ing . the effective Ha miltoni- 
ans have providedlI3llJOl22'l^l^ consid- 
erable insights into the nature of the eigenstates and IVR- 
For example, the genesis of the isomerization statedly in 
HCP and determiningjij the counter-rotation modes in- 
volving the Hydrogens in C2H2 have all been possible due 
to such detailed quantum-classical correspondence stud- 
ies. However, most of these detailed studies are restricted 
to two dimensional or quasi-two degree-of-freedom sys- 
tems with littleEJ or no work on three or higher degree of 
freedom systems. The reasons for this are many fold and 
mainly have to do with the fact that correlating quantum 
eigenstates with classical .Epmcare surface of sections, pe- 
riodic orbits, and Husim£H distribution functions from a 
phase space perspective is considerably harder for three 
or more degrees. From an eigenstate perspective, it is 
clear that identifying relevant couplings in a particular 
energy range is intimately linked to the possible dynam- 
ical assignments of the highly excited eigenstates. 

In principle, the required information is buried in the 
splitting pattern as revealed by a eigenstate resolved 
spectrum. However, in practice the task of analyzing 
the splittings in order to determine even the nature of 
the first tier states is nontrivial. In the happy instances 
when the bright state is perturbed strongly by a single 
state, there is a two-state scenario, and theoretical as 
well as experimental techniques can be brought to bear 
on the problem. Examples of such two state interac- 
tions can be found in the experiments by Boyarkin eL. 
al. on the IVR from O-H stretch overtones in methanols 
and Mork et. -al. on vibrational mode-coupling in 1,2- 
difluoroethane ^ j— .AZatious level selection schemes have 
been proposedoEJu'EII which allow one to determine the 
dark states tier structure for a specific bright state of in- 
terest. However, as soon as the two state picture gets 



superseded by a situation where two or more states are 
coupled strongly to the bright state ths-analysis is quite 
involved. A previous theoretical studyEa on the assign- 
ment of the highly excited vibrational states of H2O pro- 
vides a good example of the complications that can po- 
tentially arise even in a three-state situation. The com- 
plexity of the analysis escalates with increasing density of 
states and with the presence of nonperturbative quantum 
effects like dynamicaltij and chaos-assisted tunnelingc3. 
It is important to note that recent high .resolution ex- 
periments using double resonance methodscirHLJcJ have 
made it possible to elucidate the nature of the dark states 
that are coupled to the zeroth-order bright state. Ana- 
lyzing the double resonance spectra with the technique of 
the hierarchical tree analysisEJO, introduced by Davis, 
and the traditional perturbation methods yield impor- 
tant information about the nature of the dark states and 
the number of tiers participating in the IVR process. In 
addition, with considerable effort, it is possible to identify 
the nature of the coupling i.e., anharmonic or copplis, be- 
tween the bright state and the first tier statesuJ. Such 
an approach leads to a certain amount of insight on the 
important energy flow pathways from a bright state of in- 
terest. Nevertheless, identifying a dominant anharmonic 
or coriolis perturbation among the various possibilities is 
still an important and difficult problem. 

One of the principal aims of this work is to develop a 
systematic and experimentally accessible method to de- 
cide the relative importance of perturbations in a given 
energy range. It is also desirable to be able to an- 
alyze multidimensional effective spectroscopic Hamilto- 
nians in terms of the classical-quantum correspondence 
without explicit determination of the periodic orbits. 
In particular, this work begins to provide a theoreti- 
cal foundation for the central feature that has emerged 
from the various experimental and theoretical studies 
of IVR, i.e. the strong correlation between the eigen- 
value and eigenfunction fluctuation measures. Manifes- 
tations of such correlations are evident in the hierarchical 
tier structure, power law decay of. the survival proba- 
bility at intermediate timesaHoEj, nonejtgodic energy 

and local- 




flow even in large molecules like Benz 
ized eigenstates at faidy high energiesE3'E3'lil. The hier- 
archical tree analysisEj was one of the first attempts to 
take into account this correlation in a systematic man- 
ner. More recently— Gruebele has addressed the prob- 
lem by introducing^] the matrix fluctuation-dissipation 
(MFD) method. Apart from having numerical advan- 
tages, the MFD theorem explicitly highlights the corre- 
lation between spectral intensities and eigenvalues. Per- 
haps more relevant in the context of the present work 
is the important role played by the parametric varia- 
tion of the eigenvalues with changing bright-dark cou- 
pling strengthEj. Interestingly, many years ago Weiss- 
man and Jortner discoverectj a remarkable correlation 
between the topological features of the quantum Poincare 
maps (basically the Husimi distribution function) and the 
sensitivity of the energy levels to the strength of the non- 
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linear coupling. In what follows, it will become clear that 
the current work provides a quantitative and qualitative 
basis for the observations by Weissman and Jortner. 

Independently, considerable efforts have been made by 
the quantum chaos community towards understanding 
the origins and mec hanism s o f-e igen s ta te localization in 
fully chaotic systems^ESE3e3'E3'E3EilE3. In this extreme 
limit, quantum ergodicity is conjectured to generate, ca- 
sults consistent with random matrix theories (RMT)E3a. 
A key feature of RMT is the complete absence of corre- 
lations between eigenvalue and eigenfunction properties. 
Recently one of us (ST)E3 introduced a sensitive measure 
of eigenstate localization which is based on correlations 
between intensities and eigenlevel velocities. Eigenlevel 
velocities are more properly defined as the slopes of the 
level curves generated by changes in some system-param- 
eter. It was demonstrated in a series of papersE§E3E!l that 
this intensity-level velocity correlation function provided 
an ideal measure to explorCr-the systematic deviations 
from the predictions of RMTE3 due to localization. For 
a detailed introduction to the localization aspects and as 
to how it relates to the deviations from RMT predictions 
we refer the reader to Rcf. An attractive feature of 
a measure based on the level velocities stems from re- 
cent works demonstrating the fingerprints of a nonlinear 
resonance in the level dynamicsEJ and the ability of the 
level velocities to-distinguish between different regions of 
the phase spaceE3. In— a more general context, the cur- 
rent and recent worksEa provide a natural framework to 
analyze the spectroscopy of molecules in external fields. 
Recent examples of such "parametric spectroscopxii-jex- 
periments include the molecular Stark effect studyOE-3 on 
2-propyn-l-ol and the Laser Induced Fluorescence (LIF) 
study on NO2 in an external magnetic fieldO. These ob- 
servations provide us with a strong motivation to explore 
the potential usefulness of the correlation measure in the 
study of IVR. 

The paper is organized as follows. A brief introduc- 
tion to the intensity-level velocity correlation function 
is provided in section II. The main results for the case 
when the underlying classical dynamics is fully chaotic, 
and the RMT expectations are briefly discussed. In sec- 
tion III, the semiclassical theory for the correlation func- 
tion is constructed for classically integrable systems. Un- 
derstanding the integrable case is a necessary compo- 
nent of understanding realistic spectroscopic Hamiltoni- 
ans that often possess dynamics more closely related to 
the integrable, near-integrable, and/or mixed phase space 
regimes. The near-integrable regime is characterized 
by the introduction of resonances, and the applicability 
of first order classical perturbation theory. The mixed 
phase space regime is characterized by the co-existcncc 
of near-integrable motion and significant chaotic motion 
on the same energy surface; i.e. different initial con- 
ditions lead to different dynamical behaviors. Explicit 
expressions are given that can be verified in great de- 
tail, at least in some cases. The T% scalings of the vari- 
ous quantities are deduced as a function of numbers of 



degrees of freedom, and shown to be different from the 
chaotic case. In section IV, the intensity-level velocity 
correlation is computed for single resonance (classically 
integrable) and multiresonant (classically nonintegrable) 
two degree-of- freedom Hamiltonians. Section V provides 
a brief summary and concludes. 



II. INTENSITY-LEVEL VELOCITY 
CORRELATION FUNCTION : BACKGROUND 



Consider the effective spectroscopic Hamiltonian 
-ff(a, at;r) = H {n) + Hi (a, a* ; r) with r representing 
some parameter or parameters of interest. The eigen- 
states and eigenvalues of H will be denoted by |a(r)) and 
E a {j) respectively. In what follows, it is not necessary 
for this particular choice of the form of the Hamiltonian. 
Thus, r could equally well be a parameter of the zeroth- 
order Hamiltonian. The choice is motivated by a basic 
interest in determining the localization (or lack thereof) 
of eigenstates and hence the perturbations which have 
a large effect in a given energy range. The eigenstates 
and dynamics associated with Hq are assumed to be well 
understood which is the generic situation in IVR stud- 
ies. Throughout this article we will focus on vibrational 
dynamics and ignore the rotational effects. Nevertheless, 
the basic theory presented below can be applied to the 
full rovibrational Hamiltonian. Using the Heisenberg cor- 
respondence rule,E3 it is possible to associate a classical 
Hamiltonian H(I, <p\ r) = -ffo(I) + H\ (I, </>; r) with the 
quantum H . The form of the spectroscopic Hamiltonian 
makes it especially convenient to deal with the classi- 
cal Hamiltonian in terms of the action-angle variables 
corresponding to the zeroth-order Hamiltonian. In ad- 
dition, the classical Hamiltonian is a nonlinear resonant 
Hamiltonian representing the various resonant perturba- 
tions which couple the zeroth-order states leading to the 
flow of vibrational energy. Depending on r and the en- 
ergy of interest, the classical dynamics can exhibit the 
full dynamical range from integrable to mixed to chaotic 
behaviors. 1— in in 

Given a special state \z), the strength functior£!j|£§L3 
is defined as 



S z (E,t) 



1 



dte lEt/h {z\e- tS( - T ^/ h \z) (la) 
(lb) 



2nh 

J2p*<*(t)6(E - E a (r)) 



with p za (T) = \( z \ a ( T ))\ 2 being the overlap intensity. 
The strength function had been introduced and studied 
earlier by HellerO in the context of quantum localiza- 
tion due to eigenstate scarring. The choice of the special 
state is determined by the system of interest. From a 
spectroscopic viewpoint, the special state \z) could be 
the bright state, and then the strength function is the 
observed spectrum. Alternatively, one can choose a co- 
herent state corresponding to the phase space image of a 
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certain bright state of interest. In such a case, the inten- 
sities are essentially the Husimi function. In Sect. IV, 
we will primarily choose the latter representation in or- 
der to highlight the classical-quantum correspondence of 
the system. 

Beyond scarring, a number of other factors have been 
found to contribute towards localization effects. Ear 
instance, transport barrierSpLike broken separatricesEO, 
cantorS or diffusive motionEa in phase space have been 
linked to localization. One of the hallmarks of localiza- 
tion is the non-democratic response of the system to para- 
metric variation. Specifically the 'evolution' of the energy 
levels of the system, as measured by the level velocities, 
upon varying a parameter will not be universal. Thus, 
the independence of eigenvalue and eigenfunction fluctu- 
ation measures, a central feature of RMT, is violated. 
An analysis of the correlation between level velocities 
dE a (T)/dr and the intensities p za (T) leads to a better 
understanding of the deviations from universal behavior 
and gives insights into the nature of any localization. It 
is only very recentljCj that such a overlap intensity-level 
velocity correlation has been proposed and shown to be 
a very sensitive measure. This correlation function is de- 
fined as 

1 / s dE a (r) \ 
C z r) = — — ( Pza (r) ) (2) 

where a 2 and a\ are the local variances of the overlap 
intensities and level velocities, respectively. The brack- 
ets denote a local energy average in the neighborhood 
of E. As defined, C 2 (r) weights most the states which 
share common localization characteristics, and measures 
the tendency of these levels to move in a common di- 
rection. Within RMT, it is possible to showE3 that 
C z (t) = ± N- 1 ' 2 with N being the effective num- 
ber of states in the selected energy range. It is impor- 
tant to note that the above result for ergodic systems is 
true for every choice of \z). If there is localization, then 
there is at least one \z) for which C z {t) deviates from 
zero in a statistically significant way. The usefulness of 
C z (t) has bae*i convincingly demonstrated for the sta- 
dium billiarcEl and the baker's mapEZI. A fairly detailed 
semiclassical theory for the correlation function when the 
underlying classical dynamics is completely chaotic has 
been developedEj. That which is lacking is the theoreti- 
cal foundation for integrable, near-integrable, and mixed 
phase space systems. 

III. SEMICLASSICAL THEORY FOR C z {t) IN 
THE INTEGRABLE CASE 

The prime motivation of introducing and studying 
C z (t) for the chaotic case has to do with uncovering var- 
ious localization features of the system and hence devia- 
tions from RMT predictions. Evidently, any localization 
of eigenstates in the chaotic system is subtle and some- 
what surprising at first sight. However, localization is the 



rule rather than the exception for the integrable to mixed 
phase space regimes. Therefore, it is to be expected to 
find statistically significant nonzero values for C z (t) in 
such cases. 

It is demonstrated ahead that the correlation mea- 
sure extracts useful information for the non-fully chaotic 
regimes as well. The utility stems from the fact that 
a perturbation could non-democratically couple certain 
states or class of states leading to avoided crossings or 
like level movements within the particular class. Such a 
preferential coupling is precisely indicated by C z (t) irre- 
spective of the nature of the underlying classical dynam- 
ics. In context of the present paper, the effective spec- 
troscopic Hamiltonians typically exhibit a mixed phase 
space. As mentioned previously, the classical limit of the 
spectroscopic Hamiltonian is a nonlinear resonant Hamil- 
tonian, and the various nonlinear resonances manifest 
themselves in different energy ranges. It is quite pos- 
sible for a molecule in a certain energy range to be domi- 
nated by a single Fermi resonance. Single resonant classi- 
cal Hamiltonians are integrableLa, and. therefore one has 
conserved quantities, called polyadalj, associated with 
the motion. In a different energy range, more than one 
independent resonance can be important, and classically 
one has a nonintegrable system. The entire range of be- 
havior from integrable to chaotic is therefore possible. 
Analyzing a spectrum with C z (t) provides one possible 
route to deciphering the important resonances operative 
in the energy range of interest. Moreover, as indicated 
ahead, there are many other insights, such as indications 
of the effective numbers of degrees of freedom found in 
the ?L-scalings, that one can gain from a study of the cor- 
relation function which makes it an attractive candidate 
for exploring IVR in isolated molecules. 

It is important to note that the RMT estimate is only 
valid for the extreme fully chaotic case. In addition, given 
the richness of the quantum-classical correspondence for 
the effective spectroscopic Hamiltonians, it is both neces- 
sary and useful to come up with a semiclassical theory for 
C z (t) in such regimes. This section of the paper closely 
follows the earlier papertj concerned with the correla- 
tion function for chaotic systems. In the current work an 
analogus semiclassical theory for the integrable cases is 
provided. The theory for the ngar-integrable case will be 
provided in a later publication^. 



A. Level Velocities 

The purpose of this section is to derive explicit expres- 
sions for the mean and variance of the level velocities. 
From these expressions, the ?i-scaling properties can be 
deduced. The analogous expressions and scalings for the 
chaotic case have been derived earlierES The analysis in 
this section is important from a classical-quantum cor- 
respondence perspective as well as due to the fact that 
it can be translated into a density of states scaling. To 
focus on the integrable case we consider a single n : m 
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resonant Hamiltonian of the form: 

H(I, <{>) = H (I) + Ty/I?I% cos(m0i - n4>2) (3) 

It is well known that the above Hamiltonian is inte- 
grable because of the existence of a conserved quantity 
/ = (n/m)Ii + I2 in addition to the energy E. This 
constant of the motion I is called the polyad number. 
Since H(l, <p) is integrable one can imagine a canoni- 
cal transformation (I, <f>) — > (K, <&) in terms of which 
H = H(K). Stated differently, the variables (K, $) 
are the action-angle varibales for the single resonance 
Hamiltonian. Note that the explicit determination of 
the action-angle variables for an arbitrary single res- 
onance—Hamiltonian is very difficult if not impossible. 
JoyeuxO has given the explicit solutions for resonances 
with m + n < 4 in terms of elliptic functions. How- 
ever, not being able to construct the explicit transforma- 
tion does not preclude the existence of the action-angle 
variables for an integrable system. As pointed out by 
Joyeux, with a convenient choice of the canonical trans- 
formation, it is possible to identify one of the actions K 
with the polyad number I. Given H(K), it is possible 
to define the frequencies Vk5(K) = f2(K). The single 
resonance systems have a slight complication in the sense 
that the associated phase space has a separatrix. This 
separatrix, which separates resonant and nonresonant re- 
gions in the phase space, cannot be transformed away by 
a canonical transformation. This is reflected in the fact 
that different action variables K have to be used in the 
two regions. The separatrix structure will be ignored by 
assuming that only a local region of the phase space is 
being investigated in any given case. A more rigorous 
treatment by taking into account separatekgs would in- 
volve near-integrable semiclassical theorylZ§Ea, and that 
is left for another workllj. 

In order to determine an expression for the level ve- 
locities in terms of classical quantities and h, it is con- 
venient rte|-^faft from the smoothed spectral staircase 
functioned' 



N e {E,r) = ^f? £ [£-£ Q (r)] 



(4) 



where e is an energy smoothing term which will be taken 
to be smaller than the mean level spacing. Differentiating 
the staircase function with respect to the parameter r 
and taking the energy average yields the identity 



dN e (E,- 
d7~ 



dE a (r) 
dr 



(5) 



(Weyl) term and an oscillating term 

N e (E,r) = N{E,T) + N e (E,r) 



(7) 



For an integrable system, the oscillating part according 
to the theory of Berry and TaboiEj can be written in 
terms of the rational periodic orbits M of H(K) as: 



N € {E,t) 



2 



x exp 



£1 x — ^ 

TT)72 A M (E, t) sin[<7 M ] 

M 

eT M (E : r)' 



(8) 



It has been assumed that the curvature of the energy 
contour is approximately constant locally and in the 
sum the M = term is excluded. The rational pe- 
riodic orbits correspond to tori with rotation numbers 
otj t k = Qj/Clk — Vj/^k with [Xj,(ik being coprime inte- 
gers, fi specifies the primitive periodic orbit and M = rfi 
denotes a specific topology-of the closed orbits with r be- 
ing the repetition numberEa. In the above expression for 
N e , the amplitude and phase factors are given by 



A M (E, 



1 



27r|M|( d + 1 )/ 2 | KM | 1/2 
S m (E,t) it tt 
— g VM + A Pm 



(9a) 
(9b) 



where km is the scalar curvature of the energy contour 
H(K) = E. The action corresponding to an orbit of 
topology M is S'm = 27rK M • M, the phase r] M — M r] 
with rj being the Maslov indicies and /3m is equal to the 
sign of the determinant of the curvature matrix. The 
period of the corresponding orbit is denoted by Tm = 

The oscillatory contribution to the velocity average can 
be written down and evaluated in the h — > limit as 



dN e 



x exp 
2 



M 



fdS M 

M \ — I COS[CT M J 



V dr 



n (d+i)/ 



x exp 




i)72 J2 Am 

M 



fdS, 



M 



V dr 
(cos[(tm]) b 



(10) 



where d(E, r) is the mean level density. Similarly, as- 
suming a nondegenerate spectrum, it can be shown that 



' dN e (E, t 
. d7~ 




d(E, T ) / fdE a ( T y 



2-KE 



\ dr 



(6) 



It is possible to write down a semiclassical expression for 
the spectral staircase function as a sum of an average 



The above result follows from the argument that for a 
given orbit the actions and their parametric derivatives 
are smooth functions of energy. In the semiclassical limit, 
the derivatives are uncorrelated from the phase term 
which averages to zero in that limit. Another way to see 
this is to observe that in performing the energy averag- 
ing by stationary phase the dominant contribution comes 
from orbits satisfying cr M (i?) = which correspond to 
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the zero period orbits. By construction, the oscillatory 
part of the staircase function does not contain the zero- 
length orbits and hence the average is vanishingly small. 

Using the expression for N e and assuming the different 
orbits to be uncorrelated from each other in the semiclas- 
sical limit, one obtains 



tori. Numerical results in the next section, obtained from 
application to integrable single resonance Hamiltonians, 
confirm the derived ^-scaling for d\. The full expression 
has been verified for the rectangular billiard, but only 
the scalings will be presented in this paper. 



'dEjj) 
. dr 





where a diagonal approximation has been used for the 
sum over rational tori. The off-diagonal terms will con- 
tribute to the variance but the dominant contributions 
and ^-dependence in the semiclassical limit is expected 
to derive from the diagonal-term. For integrable systems, 
it is reasonable to assumeH the relation 




C(£,r)T 2 



(12) 



The distinction between T 2 and T (relevant to the chaotic 
case) in the above expression is analogous to the dis- 
tinction between ballistic and diffusive behaviors, respec- 
tively. In analogy to the chaotic case, we invoke the ap- 
propriate Hannay-Ozorio sum ruleE-3 



M 



h a d 
~2T 



dT 

J>2 



(13) 



to perform the sum over the rational tori. As a result we 



obtain the dimensionless variance a E as 



^2 




C(e,t)cP ~ n 



-2d 



(14) 



The relation between the mean square action changes 
and the mean square level velocities is straightforward 
involving only the density of states to complete the con- 
nection. The Ti dependence drops out easily since the 
mean level density d ~ h~ d for a e?-degree-of-freedom sys- 
tem. This is to be compared to the scalingp^ ~ ftT^ d+1 ^ 
previously established for chaotic systemsE37 Note that 
the different fi scalings are consistent with the notion 
that chaotic systems' spectra fluctuate less and are more 
rigid (constrained) than their integrable systems' coun- 
terparts. Furthermore, both the effective numbers of 
degrees of freedom and the nature of the dynamics is 
reflected in the scaling laws. This is a general feature 
arising over and over. As a consequence of the above 
result for integrable systems, the variance of the level 
velocities scales as a% ~ 0(h), and is equal to the vari- 
ance of the action changes averaged over the appropriate 



B. Overlap intensities 

In this section a semiclassical expression for the 
strength function is derived. The derivation is firmly 
rooted in the action-angle space and consistent with the 
Berry- Tabor rational periodic orbit formalism. The h 
scaling of the intensity is derived and shown to be differ- 
ent from that of its chaotic counterpart. 

As with the level velocities it is possible to decompose 
the strength function S z (E 1 t) into an average and an 
oscillatory part 



S z (E,t) = S z (E,t) + S z (E,t) 



(15) 



where \z) is a coherent state and p za = \(z\a(T))\ 2 . The 
average or the smooth part is the contribution from zero 
length trajectories. It can also be viewed as the Fourier 
transform of the extremely rapid initial decay due to the 
shortest time scale of the dynamicsO. 

The average intensities can be written in terms of the 
smooth part of the strength functions as: 



<P*«(t)) c 



1 



S z (E,t) 



(16) 



The quantity of main interest is thus 
1 



S z (E,t) = 



{2irh) d 



A W (K,$)5[E - H(K)]dKd<i> 

(17) 



A W (K, <&) is the Wigner transform of a Gaussian 
wavepacket centered at (K, <&) and can be represented 
locally in the action angle space (K, <&) by 



A W (K,*) = 2 d exp 



-(A*) 2 - i(AK) ; 



(18) 



The vectors A$ and AK have components (A3>)j = 
(<f>j - $j)/crj and (AK)j = (Tj(Kj - Kj) respectively 
where a 2 = h/2Kj. In the above expression, a dimen- 
sionless variable has been scaled to one in order to have 
equal uncertainty scalings in the angle and action vari- 
ables. As the Hamiltonian is independent of the angles, 
the angle integrals can be performed in the semiclassi- 
cal limit. A change of variables hz = AK leads to the 
expression 



S z 



1 



T d/2 



J dze- z2 5 [E-H (Tiz ■ a' 1 + K)] (19) 
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whose evaluation gives: 



dze 



T d/2 



-1/2 



3 = 1 



(20) 



The final form results from expanding the Hamiltonian 
about K and evaluating the integral at E = H(K), which 
thus pertains to the peak of the strength function. The 
S-scaling can be extracted giving 



(p»a(r)> ~ ^ 1/2 



(21) 



It is important that the derived ?i-scaling for the inten- 
sities is independent of the nature of the dynamics, i.e. 
integrable or chaotic, as must happen. 

The oscillatory part of the strength function is a dy- 
namical object and can be expressed in the action-angle 
space as 



This is consistent with the fact that we have used a local 
form of the coherent state in the action-angle coordinates. 

The Gaussian localization in action-angle space implies 
that the dominant contribution to S z arises from trajec- 
tories which are closed on the torus K i.e., $ « = $. 
In other words, the dominant contribution comes from 
rational tori Km- This motivates an expansion of the 
action as 

S K (*, *'; E) a S M (*, *; E) + K • (* - *') (26) 

In the above expansion the quadratic terms have been 
neglected. As discussed earlier, the rational torus action 
Sm = 27tKm • M. This is reasonable for an integrable 
system since the initial Gaussian in phase space simply 
shears with time. Using the above expansion of the ac- 
tion, the integrals over the angles can be performed by 
stationary phase giving 

~ Sz{E ' T) * ft (2 t"i-i)/^ |I3M|1/acos[gM] 



M 



xe 



-(AK M ) 2 /t 2 



(27) 



SJE,t) = —ha I d*d*'(z|*)G($,*';£;)(*'|z) 

7T 



(22) 



The above result for the oscillatory part of the strength 
function can be used to obtain the variance of intensities 



where G(«&, <&'; E) is the energy Green's function—The 
semiclassical expression for the Green's function isL3 



1 V^ln 1 1/2 

ih(2Trih)( d - 1 )/ 2 *g K| 



xexp( -5 K (*,*';S) 



(23) 



where the sum is over all possible trajectories (associated 
with tori K) that go from to 3> at an energy E and 
their repetitions. Sje(&, E) is the action associated 
with a trajectory (including the repetitions) and Z?k is 
the usual stability determinant of order (2d + 1) x (2d+ 1) 
given by 



lire ~ 

-^(S 2 z (E,r)) E (28a) 

\ M 



d(2ir) d - 2 h 

xe -2eT M /h\ 
1 1 



(28b) 



where we have used the diagonal approximation for the 
sum over M. Since the trace of the Green function is 
the density of states given by Berry and TaborEj and the 
density of states is the energy derivative of the spectral 
staircase function, then |Z?m| = Tmp^m/ '(27r) d_1 . Hence, 
given the Hannay-Ozorio sum rulelil for £ M A Ia stated 
in the previous section, here 



D K = 



0B9#' dE 2 



(24) 



E 

M 



D 



M 



h d d 

{2n) d 



dT 



(29) 



The product of the wavefunctions in the coherent stata 
basis can be obtained via an inverse Wigner transformEj 
of the phase space Gaussian. The inverse transform is 
evaluated by stationary phase to obtain 



<*I*X*'I*> 



1 



d/2 rr exp 

Ik •(*-*') 



-i(A*) 2 -i(A*') 2 



(25) 



It follows that 



n, - 



2 

3 "3 



Jld-1 



-2(AK M ) 2 /t : 



M 



(30) 



which again has been verified for the rectangular billiard, 
but will not be shown here. The leading order ?i-scaling 
of the intensity variance is 



2 %a 
tz ~ n 



(31) 



The evaluation by stationary phase is valid if the location 
of the Gaussian is away from the edges of the phase space. 



In the chaotic caseE3, the RMT expectation is a 2 ~ Ti 2d 1 
which follows from the scaling of (p za )- The intensities 
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follow a Porter-Thomas distribution, which has only one 
parameter given by its mean; the variance is always twice 
the square of the mean. Again, the integrable limit scal- 
ing of the variance of the intensities is different from the 
chaotic limit which is a manifestation of the differing na- 
ture of the underlying classical dynamics. It follows that 
because the variance is on a much greater scale, the dis- 
tribution of strengths is necessarily singular for the inte- 
grable case in the limit of Ti —> (i.e. mostly near zero 
values except for a few extremely large values). 



C. Correlation function 

The semiclassical theory for overlap intensities and 
level velocities outlined in the previous subsections can 
now be used to construct the correlation function. Be- 
ginning with the covariance, it can be shown that 



Cov z (r) 



Pza(T) 



dE a (r) 
dr 



(32) 



The semiclassical expressions for the strength function 
and the parametric derivative of the staircase function 
give 



/2 



OS 



M 



dr 



3 -(AK M ) 2 /fi 2 e -2eT M /B 



Again, the diagonal approximation has been made in the 
above expression. Following arguments parallel to those 
in the previous section leads to the expression 



x e -(AK M ) ! /6 2 \ 

/ M 



(34) 



where th = h/e is the Heisenberg time. In order to 
make further progress in understanding the correlation 
function, it is necessary to know the time dependence of 
the M-averaged expression. It is important to note that 
the Gaussian weighting factor will not decouple from the 
paramteric action derivative at long times. Assuming 
that the parametric derivative <9Km/<9t to be a weak 
function of the period gives the approximation 



dS 



M 



M 



leading to the final result 



Cov z (t) 



T d~l 



(35) 



(36) 



Combining this with the results of the previous subsec- 
tions generates the semiclassical expression for the cor- 
relation function. It is 



C,(r) 



C(r) 



CV2(£,r)(e- 



2(AK M )Vft- 



1/2 
M 



(37) 



Interestingly, there is classical information represented 
in the correlation function that is not in the covariance; 
therefore considering the separate "components" of the 
correlation function (or considering both the covariance 
and correlation functions) gives the greatest information 
about the system. In the semiclassical limit, the above 
result implies that Cov z (t) ~ h d l 2 for integrable sys- 
tems whereas the correlation C z (t) ~ h° to leading or- 
der. This last result had to emerge from the semiclassi- 
cal calculations. Otherwise strong localization could not 
persist in integrable systems in the limit as h — > 0, and 
multidimensional WKB theory could not work. 



IV. APPLICATION TO A MODEL EFFECTIVE 
HAMILTONIAN 



In order to illustrate the usefulness of the intensity- 
level velocity correlation measure introduced in the previ- 
ous sections, consider a simple two dimensional effective 
Hamiltonian 

H(n,aL,sJ) = ir (n) +T 11 5 1:1 (a ! a t ) + r 21 # 2: i(a, a f ) 



(33) 



(38) 



where -ffo(n) is the zeroth-order, anharmonic, vibrational 
Hamiltonian given by 



H (n) = UiNi + uj 2 N 2 + xnNl + a^iVf 



(39) 



The zeroth-order frequencies of the two modes are de- 
noted by uji , u>2- The mode occupation numbers are 
Nj = (rij + 1/2) and the mode anharmonicities are de- 
noted by Xjj. The full Hamiltonian is multiresonant in 
that the zeroth-order anharmonic modes are coupled by 
a 1 : 1 and a 2 : 1 resonant perturbation defined in terms 
of the operators a, at as 

i?i : i(a, a*) = a\a 2 + alai 

flaa (a, at) = a\a 2 a 2 + a\a\a\ (40) 

The parameters th and t 2 \ measure the strength of 
the corresponding resonant perturbations. The above 
spectroscopic Hamiltonian H arises in many instances 
in molecular spectroscopy. For example, with t 2 \ = 
and fop-,identical oscillators, the effective Hamiltonian 
modelsc3 the stretching modes of H2O. On the other 
hand with t\\ — 0, the effective Hamiltonian models the 
Fermi resonant interactions between the stretch and the 
bend modes in HCX3 moleculealH. We have chosen a 
multiresonant Hamiltonian to model systems where two 
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zeroth-order modes of interest can get coupled by one or 
both of the resonances depending on the pertinent en- 
ergy range. The parameter values chosen for this work 
are u>i = 1.0, lu 2 = 0.8, i u = -0.04, and x 2 2 = -0.02. 

The classical limit Hamiltonian corresponding to the 
quantum H can be easilw-,constructed using the Heisen- 
berg correspondence rule£3 and shown to be the following 
nonlinear, multircsonant Hamiltonian: 

H(1, 0) = (Jih + ijJ2h + xnlf + X22I2 
+ 2r n (hh) 1/2 cos(fa -fa) 
+ 2 T2 i(hli) 1/2 cos(fa-2fa) (41) 

It is well known that with either tu or T21 equal to zero 
the classical Hamiltonian is still integrable whereas with 
both the coupling strengths nonzero, the system is non- 
integrable. For the choice of the parameter values made 
above and noting that the model Hamiltonian is two di- 
mensional, the density of states is expectedly small. In 
order to simulate a high density aLstates situation we will 
make use of the following scalingc3 to reduce the effective 
value of Ti — ► h/c: 

LOj > LUj LOj J C, Xij ► Xij — Xij J C , 

Tn -> rn = nx/c, T 2 i -> T21 = r 2 i/c 3/2 (42) 

Note that for most of the results presented in this work we 
have used U = 1. Elsewhere, the scaling will be explicitly 
given. 

In computing the correlation functions shown below, 
the level— .velocities have been constrained to be zero- 
centerecE3. This is achieved by subtracting out the mean 
level velocity, in the energy range of interest, from the 
individual level velocities; this results in no loss of gen- 
erality. Also, the energy window is chosen such that 
the classical dynamics are essentially the same through- 
out the range. This is. essential since any bifurcation in 
the system will alteillj the level velocity structure. The 
problem of qualitative changes in the correlations due to 
the various bifurcations will not be addressed here. The 
"bright state" \z) is chosen to be a coherent state cen- 
tered on (J, -0) with J = (Ii — h)/2 and ip = (fa — fa). 
This choice of \z) is convenient for studying the system's 
classical-quantum correspondence via comparisons to the 
appropriate classical Poincare surface of sections. The 
details for the construction of the surface of sections and 
deteipiining the intensities p za can be found in an earlier 
workEJ. The level velocities are computed by invoking 
the Hellman-Feynman theorem and given by 



Note the important point that the level velocities them- 
selves are basis independent objects. The energy range 
used for computing the correlation is fixed to be [8,9]. 
For this choice of the energy range, the eigenvalues and 
eigenstates are well converged in the zeroth-order number 



basis. The choice of the energy range is also motivated by 
the fact that the two resonances overlap without leading 
to large scale stochasticity in the classical phase space. 



A. Example of integrable cases 

The multiresonant Hamiltonian has two integrable, 
single resonance subsystems that can be obtained by set- 
ting m 7^ 0, T21 = and t\\ — 0, T21 7^ 0. In the former 
case one obtains a 1 : 1 system and in the latter case 
a 2 : 1 system is obtained. Studying these integrable 
subsystems serves two main purposes. First, the evalu- 
ation of the correlation function C z (t) helps provide a 
baseline for the later interpretations in the nonintegrable 
regime. Secondly, the ?i-scalings predicted for a a and <je 
can be verified. In principle, the full semiclassical expres- 
sions could be verified. However, that would require the 
calculation of a large set of periodic orbits which would 
be time consuming for this system; the expressions have 
been verified for the simpler system of a rectangular bil- 
liard. Note that at the outset it is expected that the 
correlation will be significant for many choices of the co- 
herent states \z). Our concern, however, is to determine 
the nature of C z (t) for the integrable systems and com- 
pare it to the classical phase space structures. 

Throughout this section the nonzero coupling strength 
is set equal to 10~ 2 . In Fig. ([!]), we show the classical 
(J, ip) surface of section at an energy E = 8.5 for the 
integrable 1 : 1 and the 2 : 1 cases, respectively. The ap- 
propriate resonant and nonresonant regions for the two 
cases can be clearly identified. In the energy range of 
interest, there are no bifurcations in the classical phase 
space. In Fig. (||), we show the intensity and dimension- 
less level velocity variances as a function of h and h~ 2 , 
respectively (so that the points lie on a straight line), 
for the two subsystems. In the 1 : 1 case, the coherent 
state is located at \za) = (0, —1) in the (J, tp) surface of 
section whereas in the 2 : 1 case, the coherent state is lo- 
cated at \zb) = (—0.5,-3.75). The particular choices of 
the coherent states are made such that they are located 
inside the resonance zones and close to the appropriate 
stable periodic orbits. The agreement with the theoret- 
ical predictions is fairly good and other choices of \z) 
yield similar results except near the edges of the (J, tp) 
plane. Although not shown here, the scaling of the aver- 
age intensities (p za ) is also consistent with the theoretical 
values derived in the previous section. 

We now turn our attention to the intensity-level ve- 
locity correlations and compute them as a function of 
the coherent states spanning the surface of section; each 
correlation result is plotted at the phase space center of 
the coherent state being used, and which thereby gives 
a full phase space portrait of the correlation function's 
behavior. In Fig. (g), we show the positive and negative 
parts of the correlation C z (rn), respectively. The prime 
reason for showing the positive and negative contribu- 
tions of C z (th) separately is that this directly reflects 
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FIG. 1: Poincare surface of sections for the integrable subsys- 
tems at energy E — 8.5. The top panel corresponds to the 
1 : 1 case and the bottom panel to the 2 : 1 case. 



the phase space localization features of states with pos- 
itive and negative level velocities. As mentioned in the 
introduction, the relative sign of the level velocities in 
integrable and mixed systems is intimately connected to 
the localization of the states about stable and unstable 
structures in the phase space. For h = 1, the number of 
states in the energy range is 27 for both the 1 : 1 and the 
2 : 1 systems. Despite the low density of states, Fig. (||a) 
shows large positive, statistically significant values for 
C z (tii) around the stable periodic orbits. Similarly from 
Fig. (||b) , notice the large values for the negative correla- 
tion around the unstable periodic orbit; the positive and 
negative peak values are 0.92 and 0.42 respectively. The 
results for the integrable 2 : 1 case are shown in Fig. 
Once again the correlation C z {j2i) takes on large posi- 
tive (0.83) and negative (0.65) values at the stable and 
unstable periodic orbits respectively. The strong posi- 
tive peaks in the correlations are reminiscent of similar 
structures arising-, due to the bouncing ball states in the 
stadium systemE3. Although the stadium is completely 
chaotic, the analogy to the present integrable case is ap- 



qL^U 1 1 1 , 1 1 |q 

0.2 0.4 0.6 0.8 1 
h 

FIG. 2: 7i-scaling for the intensity variance (o a ) shown as solid 
circles and the dimensionless velocity variance (3\e) shown as 
open circles. The lines correspond to the theoretical esti- 
mates, (a) integrable 1 : 1 with m = 10~ 2 ,T2i = 0. (b) 
integrable 2 : 1 with rii = 0, T21 = 10~ 2 . 



propciate. This has to do with the fact that a previous 
workl23 has shown that the states localized in the center of 
the resonance-zone exhibit linear parametric motion. It 
is well knownEJ that the bouncing ball states also exhibit 
linear parametric motion apd, indeed, a WKB quantiza- 
tion of such states leads toE3 oe ~ fi~ 2 which compares 
favorably with our integrable system theory. Strong cor- 
relation at the classical fixed points have been observed 
earlieiO in the context of the bakers and the standard 
map. The present results are consistent with the previ- 
ous work, and a significant observation is that the corre- 
lation is a faithful representation of the relevant resonant 
portion of the phase space. In other words, depending on 
the parameter being varied, the corresponding region of 
the classical phase space is highlighted. A particularly 
interesting feature of the correlations in Figs. (0,0) has 
to do with a noticeable dip at the classical fixed points. 
This seems to be a fairly general feature associated with 
fixed points in the classical phase space. In a sense, the 




11 




FIG. 3: (a) Positive part of the correlation C z (rn) for the 
integrable 1 : 1 subsystem with each point representing a 
possible coherent state \z) in the (J,ip) plane. The axes are 
(J, ip) with identical limits as in the surface of section in figure 
1. The z-axis ranges from to 1. Note that all the figures 
illustrating the correlation function have the same range for 
the axes as in this figure. (b)Negative part of the correlation 
Cz(th) for the integrable 1 : 1 subsystem. 



correlation, constructed with quantum eigenvalues and 
eigenstates, is capable of locating the periodic orbits of 
the system. We do not intend to give a detailed account 
of this feature except to mention that it arises from the 
large intensity variance a a at the fixed points. 



B. Example of mixed phase space case 

Consider the full multiresonant Hamiltonian with 
Tn = 7"2i = 10~ 2 . The theory developed in this work 
is not strictly applicable in this regime, although if prior 
experience can be used as a guide, it will be relevant to 
the regular phase space regions or KAM islands (as will 
be seen ahead) . On the other hand, a mixed phase space 
is generic to many systems, and it is essential to under- 
stand the structure of the correlation function in such 
cases. It is important to note that a theoretical basis 



for the mixed phase space case is exceedingly difficuJt, 
and currently the next best approach is to attemptE3 a 
description for the near- integrable systems. 

Four surfaces of section in the energy range [8,9] are 
shown in Fig. (|5|). The classical dynamics does not vary 
significantly in this energy range and exhibits a mixed 
phase space. The two resonances have overlapped, gener- 
ating stochastic regions in the phase space. Yet, it is still 
possible to identify the stable fixed points correspond- 
ing to the two resonances. Figure (^|) shows the posi- 
tive and negative components of the correlation function 
C z (tii). As compared to the integrable case in Fig. (||), 
there are many common features including the large pos- 
itive values at the stable fixed point and the associated 
dip. This indicates that the linear parametric motion 
persists at these energies and coupling strengths. There 
are, however, important differences between the mixed 
phase space and the integrable cases. In particular the 
negative regions of the correlation in Fig. (|^b) show a 
lot more structure. The dip associated with the unsta- 
ble fixed point is hardly noticeable. This is consistent 
with the destruction of the separatrix due to the overlap 
of the two primary resonances. Notice that the nega- 
tive values of the correlation are mainly associated with 
the chaotic regions of the phase space as opposed to the 
positive values shown in Fig. (^|a). The nonuniformity 
of the negative C z (rn) implies that there exist eigen- 
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states in the energy range which are affected strongiy by 
the 1 : 1 perturbation and which share common local- 
ization features. Figure (^|a) clearly indicates the exis- 
tence of strongly localized states close to the center of 
the 1 : 1 resonance channel. Roughly speaking, although 
in the limit of strongly chaotic dynamics, the chaotic re- 
gion should show no structure, the chaotic region in a 
mixed phase space system has a great deal of structure 
related to cantori or broken separatrices, and "stickiness" 
whereby trajectories become trapped near tiny regular 
regions. Transport is not rapid and uniform, and the 
localizing effects on the eigenstates show up in the corre- 
lation function. In order for more details of this structure 
to be visible in Fig. (^|a) , a much higher density of states 
(smaller value of ft) would have to be used. 

Figure (pi) shows the analogus results for the correla- 
tion C 2 (t2i). The results are similar to that of the 1 : 1 
case. Parametric variations in T21 highlight the 2 : 1 re- 
gion of the phase space. Note that it is quite possible 
to choose a state \z) such that both the correlations are 
negligibly small. This implies that the chosen \z) is not 
strongly affected by either of the resonances. Thus, if a 
coherent state \z) lies in a region that is strongly affected 
by a certain perturbation, then the associated correlation 
will highlight the relevant part of the phase space. 

The model system considered in this paper has two 
degrees of freedom. Consequently it is possible to ob- 
tain detailed information about the phase space struc- 
ture and localization features of eigenstates. In three or 
more degrees of freedom systems the standard approach 



FIG. 6: (a) Positive and (b) negative part of the correlation 
C*(tu) for the nonintegrable system. 



of surface of sections and phase space distribution func- 
tions cannot be employed for well known reasonsL3. The 
results in this section demonstrate that the correlation 
function highlights the various classical structures (peri- 
odic orbits, tori, etc.) in the phase space. Nevertheless, 
the correlation function is not subject to the usual con- 
straints that have to be satisfied by a surface of section. 
This suggests that it is possible to compute the correla- 
tion over various reduced dimensionality spaces and gain 
insights into the system. For instance, the correlation 
could be averaged over the angle variables, and obtain 
an average correlation as a function of the actions alone. 
As a simple example in Fig. (JsJ) , we show the angle aver- 
aged correlations {C z (T))g corresponding to the results in 
Fig. (j(],|7]) . The peaks of the averaged correlation function 
are very close to the center of the respective resonance 
zones. Hence, (C z (r))g is providing information on the 
local resonance_structure associated with a particular \z) 
which is linkedEJ experimentally to the gateway states. 

As metioned in the introduction, experimentally there 
is access to very few bright states in a particular en- 
ergy range. One of the important questions, given the 
bright state and its associated spectrum, concerns the 
relative effect of the various perturbations in a certain 
energy range. It is of particular interest to demonstrate 
how the above question can be addressed with the over- 
lap intensity-level velocity correlation function. In or- 
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(a) 




(b) 



FIG. 7: (a) Positive and (b) negative part of the correlation 
Cz(t2i) for the nonintegrable system. 



der to determine the usefulness of the correlation in this 
regard, we select three different bright states \za), \zb) 
and, \zc)- The first two have been introduced in the pre- 
vious subsection, and the third one \zq) is centered at 
(J, ijj) = (0.6,-1.5) corresponding to the chaotic region 
of the phase space. The spectra associated with the three 
states are shown in Fig. In order to have a higher 
density of states, the value of h = 1/3. For the chosen 
value of h, there are 250 states in the energy range shown 
in Fig. ([}]). The correlations C z (th) and C z (r 2 i) associ- 
ated with each of the spectra are also shown. It is clear 
that for the bright state \za) the 1 : 1 perturbation has 
a much larger effect as compared to the 2 : 1 perturba- 
tion. The roles are reversed in the case of the bright state 
\zb)- In the case of the state \zc), both of the correla- 
tions are negative and small. However, the correlations 
are not insignificant since, as a rough estimate, a RMT- 
like result would give C z ~ ± 0.063. Moreover, the 
correlation is almost two times larger for the 1 : 1 per- 
turbation as compared to the 2 : 1 perturbation. ^From a 
time-dependent point of view, recurrences in the survival 
amplitude (z\z{t)) are connected to classical structures 
in the phase space. The correlation function is providing 
clues as to the precise nature of the classical phase space 
structures which are dictating the dynamics of the bright 
state of interest. 

Finally, a closer look at the correlations shown in 



U 



0.2 
0.1 

CD 

N 



-0.1 



-0.2 



/ V 


1 1 1 1 1 - 
\ 


/ A 


\ 


/ ' \ 


\ 


- / / \ 


\ 


/ / \ 


\ 


i ^ \ 


\ 


/ / \ 


\ f N. 


~ / /' 


\ \ / 


/ / 


\ V 


■ / / 


J y' 


\ r 








\ /' 






-1,1 





-3 




J 



FIG. 8: Angle averaged correlations in the nonintegrable case. 
The solid line corresponds to the 2 : 1 case and the dot-dashed 
line corresponds to the 1 : 1 correlation. The maximum in 
each case is very close to the corresponding resonant center J 
values. For the parameters chosen in this work JJi ~ —3.75 
while Jj :1 is consistent with the maximum of the averaged 
correlation. 
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FIG. 9: Spectra associated with three different coherent states 
and the correlation functions, h — 1/3 in this figure. (A) \za) 
centered close to the 1 : 1 stable fixed point. (B) \zb) centered 
close to the 2 : 1 stable fixed point. (C) \zc) centered in the 
chaotic region. The bottom right panel shows the correlations 
for the three coherent states with circles denoting C z (rn) and 
squares denoting C z (t2i). The lines are drawn as a guide to 
the eye. 
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Figs. (|| 0), compared to their integrable counterparts, 
indicates the spreading of the correlation due to the de- 
struction of barriers in the phase space. The classical sur- 
faces of sections shown in Fig. (||) have been computed 
up to times significantly larger than the heisenberg time 
th ■ The fact that the correlations are indicating "ampli- 
tude transport" through broken KAM tori or separatrices 
gives a unique probe of the transport. Moreover, the un- 
even spreading of the correlations suggests the existence 
of partial barriers. This observation raises the possibil- 
ity of understanding transport phenomena in the phase 
space via the correlation function. Further work is war- 
ranted to establish what information is contained, and 
how to extract it from the correlation behaviors in the 
chaotic regions of mixed phase space systems. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have explored the possibility of study- 
ing IVR through the overlap intensity-level velocity cor- 
relator. A semiclassical theory for the correlator has 
been developed in the case that the underlying clas- 
sical dynamics is integrable. Combined with the ear- 
lier world£3 providing the semiclassical theory for chaotic 
systems, the correlator is an ideal candidate for a de- 
tailed study of the classical-quantum correspondence for 
a general system. Application to an effective spectro- 
scopic Hamiltonian clearly supports our point of view. 
Although the model Hamiltonian used in the present 
study is two dimensional, it is evident that the corre- 
lator is well suited for multidimensional systems with a 
high density of states. It was demonstrated that choos- 
ing different parameters, corresponding to the different 
perturbations, highlighted the appropriate parts of the 
phase space displaying prominent localization features. 
The fact that these localization features are often con- 
nected to classical structures like tori, resonance zones, 
and periodic orbits provides an avenue to understand 
the relative importance of the various perturbations in 
a given energy range. From the IVR perspective this 
information is intimately connected to the dynamics of 
an initial bright state of interest. If the dynamics of a 
particular bright state is influenced strongly by specific 
classical structures, then the correlation will invariably 
reveal the relevant terms in the molecular Hamiltonian 
which play an important role. 

The manner in which we have applied the technique 
in this work is useful in situations where one has an ef- 
fective spectroscopic Hamiltonian at hand. For a specific 
molecule, the coupling strengths are fixed and experi- 
mentally one cannot vary the perturbation strengths in 
a direct fashion. However, the standard practice of iso- 



topic labeling/substitution achieves this in an, indirect, 
albeit uncontrollable, way since such studiesEJ do vary 
the nature of the dark states relative to the bright state. 
A more direct possibility would be to use various ex- 
ternal fields or parameters to generate the level motions; 
but then the connection between the fields and the possi- 
ble Hamiltonian terms would need to be explored. From 
a theoretical standpoint, the computation of the inten- 
sities requires the eigenstates which can be a daunting 
task for molecules with a large number of coupled modes. 
Significant advances in computational methods in recent 
years is certainly encouragingu. However, under some 
circumstances, a hybrid approach of correlating the ex- 
perimental intensities with the computed level velocities 
could be .tried. Another possibility is to adapt the MFD 
approached where intensities are related to the paramet- 
ric derivatives of the eigenvalues. 

Throughout this paper we have assumed the availabil- 
ity of an eigenstate resolved spectrum for the system. 
At a coarser level, corresponding to shorter time scales, 
there is the dispersed fluorescence spectrum which con- 
tains key information regarding IVR and approximately 
conserved quantities. It would be interesting to have an 
explicit time-dependent view of the IVR process i.e., de- 
ciding the perturbations which dictate the dynamics of a 
bright state over different time scales. For this purpose, 
it is appropriate ta,work with the finite time or smoothed 
strength functionO defined as: 

SP(E,t) = -L [ T dte lEt / h (z\e~* S ^ h \z) (44) 
2irh J _ T 

There is a problem with using the finite time strength 
function since the concept of eigenstates and hence that 
of level velocities is not meaningful. We are currently 
exploring the correct finite time analog of the correlator. 

As stated earlier there are many aspects of the 
intensity-level velocity correlator which warrant further 
studies. Foremost among them are phase space trans- 
port, spectroscopy in external fields and quantum versus 
classical routes to localization. Efforts are currently un- 
derway to systematically approach and understand these 
features in molecular systems. 
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